Guided Modes in Negative Refractive Index Waveguides 



O 

o 

(N 

> 
O 

in 



43 ■ 

Oh! 
i 

C/3 
C/2 ' 

<3 : 

C/3 ■ 

O ■ 

• i-H . 

C/3 . 
>v 
43 ' 



> 
in 

o 



(N 
O 

o 

>>: 

43 



X 



Ilya V. Shadrivov, Andrey A. Sukhorukov, and Yuri S. Kivshar 
Nonlinear Physics Group, Research School of Physical Sciences and Engineering, 
Australian National University, Canberra ACT 0200, Australia* 

We study linear guided waves propagating in a slab waveguide made of a negative-refraction- 
index material, the so-called left-handed waveguide. We reveal that the guided waves in left-handed 
waveguides possess a number of peculiar properties, such as the absence of the fundamental modes, 
mode double degeneracy, and sign-varying energy flux. In particular, we predict the existence of 
novel types of guided waves with a dipole-vortex structure of the Pointing vector. 
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Recent experimental demonstration of novel compos- 
ite materials with a negative index of refraction [1] opens 
up a unique opportunity to design novel types of de- 
vices where electromagnetic waves are governed in a non- 
conventional way. The history of these materials begins 
with the paper by Veselago [2], who studied the wave 
propagation in a hypothetical material with simultane- 
ously negative dielectric permittivity e and magnetic per- 
meability /i. Such media are usually termed left-handed 
materials, since the electric and magnetic fields form a 
left set of vectors with the wave vector. Already back in 
1968, Veselago predicted a number of remarkable prop- 
erties of waves in left-handed materials, such as negative 
refraction, reversed Doppler and Vavilov-Cherenkov ef- 
fects. However, structures with both negative e and fi 
have not been known until recently, although materials 
with negative dielectric permittivity are known (e.g. a 
metal below the plasma frequency). 

The study of microstructured metallic materials for 
magnetic resonance imaging [3] has shown that such 
structures can be described in terms of effective mag- 
netic permeability which becomes negative in the vicin- 
ity of a resonance frequency. It was expected that mixing 
the composite materials with negative magnetic perme- 
ability [3] with those possessing negative dielectric per- 
mittivity [4] would allow to create a novel type of meta- 
materials with a negative index of refraction. Indeed, 
both numerical simulations [5] and experimental results 
[1, 6] confirmed that such left-handed (or negative-index 
refraction) materials can be fabricated and analyzed. 

One of the first applications of the negative-refraction 
materials was suggested by Pendry [7], who demon- 
strated that a slab of a lossless negative-refraction ma- 
terial can provide a perfect image of a point source. Al- 
though the perfect image is a result of an ideal theoreti- 
cal model used in Ref. [7] , the resolution limit was shown 
to be independent on the wavelength of electromagnetic 
waves (but can be determined by other factors such as 
loss, spatial dispersion, etc.), and it can be better than 
the resolution of a conventional lens [8] . 

The improved resolution of a slab of the negative- 
refraction material can be explained by the excitation 
of surface waves at both interfaces of the slab. There- 
fore, it is important to study the properties of surface 



waves at the interfaces between the negative-refraction 
and conventional materials. So far, the frequency disper- 
sion of surface waves at a single interface and in a slab 
waveguide was calculated numerically only for particular 
medium parameters [9]. 

In this Letter, we study, for the first time to our knowl- 
edge, the structure and basic properties of electromag- 
netic waves guided by a left-handed waveguide. In order 
to emphasize the unusual and somewhat exotic properties 
of such waves, we compare them with the guided waves 
of conventional planar dielectric waveguides. We reveal 
that the guided modes in left-handed waveguides differ 
dramatically from conventional guided waves, and they 
possess a number of unusual properties, including the 
absence of the fundamental modes, double degeneracy of 
the modes, the sign-varying energy flux, etc. In partic- 
ular, we predict the existence of novel types of guided 
waves with a dipole-vortex structure of the energy flux 
and the corresponding Pointing vector. 

We consider a symmetric slab waveguide in a conven- 
tional planar geometry [see, e.g., the top left insert in 
Fig. 1(a)]. In the general slab of the thickness 2L 

is made of a material with dielectric permittivity €2 and 
magnetic permeability /12, which both can be negative 
or positive. We assume that the surrounding medium is 
right-handed, and is therefore characterized by both pos- 
itive ei and It is well known that a slab waveguide 
made of a conventional (right-handed) dielectric material 
with e2 > and \ii > creates a non-leaky waveguide 
for electromagnetic waves, provided the refractive index 
of a slab is higher than that of the surrounding dielectric 
medium (cladding), i.e. £2/^2 > However, in the 

following we demonstrate that this simple criterion can 
not be applied to the waveguides made of a left-handed 
material. 

To be specific, below we describe the properties of the 
TE guided modes in which the electric field E is polarized 
along the y axis. A similar analysis can be carried out 
for the TM modes, and these results will be presented 
elsewhere. From the Maxwell's equations, it follows that 
stationary TE modes can be described by the following 
equation for the scalar electric field E = E y , 
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E = 0, (1) 
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FIG. 1: Comparison between the conventional (a) and left- 
handed (b) guided modes of a slab waveguide. The dashed 
and solid curves correspond to the left- and right-hand sides 
of the dispersion relations in Eqs. (4) and (5), respectively. 
Intersections of these curves indicate the existence of guided 
modes. Three dashed lines in each plot correspond to waveg- 
uides with different parameters pi > pi > P3, but the fixed 
ratio /U2//U1. Inserts show the waveguide geometry, and the 
transverse profiles of the guided modes. 
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FIG. 2: Frequency dispersion curves for the three lowest-order 
guided modes of the left-handed slab waveguide (L — 2 cm). 
Insets show the characteristic mode profiles. 



where oj is the angular frequency of the monochromatic 
waves. 

The guided modes can be found as special stationary 
solutions of Eq. (1) having the following form, 



E(x, z) = E (x)e 



Ihz 



(2) 



where real h is the wave propagation constant and Eq (x) 
is the spatially localized transverse profile of the mode. 
After substituting Eq. (2) into Eq. (1), we obtain an 
eigenvalue problem that possesses spatially localized so- 
lutions for 



because only in this case the mode amplitude decays 
away from the waveguide, Eq(x) ~ exp(— |x|«i). 

We solve the eigenvalue problem in each of the ho- 
mogeneous layers, and then employ the corresponding 
boundary conditions following from Eq. (1). As a result, 
we obtain the dispersion relations which define a set of 
allowed eigenvalues h, 



(/siL) = ± — {k 2 L) tan ±1 (k 2 L), 



(4) 



where the signs (+) and (— ) correspond to the symmet- 
ric and antisymmetric guided modes, respectively, and 
k 2 = \{u 2 j c 2 )e 2 p, 2 — h 2 ] 1 / 2 . When k 2 is real, the corre- 
sponding modes can be identified as "fast waves", since 
their phase velocity uj /his larger than the phase velocity 
in an homogeneous medium with e 2 and p\, 2 . 

The parameter k 2 becomes purely imaginary for "slow 
waves" , when the propagation constant h exceeds a crit- 
ical value. Then, it is convenient to present Eq. (4) in an 
equivalent form using n 2 = ik 2 , 



{k\L) — — — (k 2 L) tanh ±1 (K 2 L). 



(5) 



Following a standard textbook analysis (see, e.g. 
Ref. [10]), we consider the parameter plane (k 2 L,KiL), 
and also extend it by including the imaginary values of 
k 2 using the auxiliary parameter plane (k 2 L, KiL). In 
Figs. l(a,b), we plot the dependencies described by the 
left-hand (dashed) and right-hand (solid) sides of Eqs. (4) 
and (5), using a relation between the parameters, 



( Kl Lf + (k 2 L) 2 



L 2 (ui 2 /c 2 ) (e 2 /j, 2 - ei/ii) = p. (6) 



ei/Lti > 0, 



(3) 



In Figs. l(a,b), we draw three dashed lines corresponding 
to different slab waveguides, having the same ratio p\j p. 2 . 
The intersections of a dashed line with solid curves indi- 
cate the existence of various guided modes. We present 
results for a conventional (right-handed) waveguide in 
Fig. 1(a), in order to compare them directly with the cor- 
responding dependencies for a left-handed slab waveguide 
in Fig. 1(b). 

First of all, the analysis of Eqs. (4) and (5) confirms 
the well-known textbook results that a right-handed slab 
waveguide can only support "fast" guided modes, which 
exist when the waveguide core has a higher refractive in- 
dex than its cladding, i.e. for e 2 fx 2 > e\ In this case, 
there always exists a fundamental guided mode, which 
profile does not contain zeros. The conventional waveg- 
uide can also support higher-order modes, their number 
is limited by the value 2p 1 / 2 /7r. These various regimes 
are illustrated in Fig. 1(a) with different dashed lines. 

However, we find that the properties of the left-handed 
slab waveguides are highly nontrivial. First, such waveg- 
uides can support "slow" modes, and they are either sym- 
metric (node- less) or antisymmetric (one zero). Such so- 
lutions represent in-phase or out-of-phase bound states 
of surface modes, localized at the two interfaces between 
right and left media. In the conventional case of both 
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FIG. 3: Surface waves in a slab waveguide for the case 
£2/^2 < tipi an d P'l < pi- Shown are (a) the propagation 
constant and (b) the energy flux vs. the slab thickness pa- 
rameter L. Solid and dashed lines correspond to strongly and 
weakly localized modes, respectively. Insets show the charac- 
teristic mode profiles. 



positive e and p, such surface waves do not exist, however 
they appear when the magnetic permeability changes its 
sign (for the TE polarization). Thus, the guided modes 
can be supported by both low-index and high-index left- 
handed slab waveguides. 

Second, the conventional hierarchy of "fast" modes 
is removed. Specifically, (i) the fundamental node-less 
mode does not exist at all, (ii) the first-order mode ex- 
ists only in a particular range of the parameter values 
p, and it always disappears in wide waveguides, when p 
exceeds a critical value, and (iii) two modes having the 
same number of nodes can co-exist in the same waveg- 
uide. We illustrate some of these nontrivial features in 
Fig. 1(b). 

Frequency dispersion of the guided waves in the left- 
handed waveguides should be studied by taking into ac- 
count the dispersion of both e 2 and p 2 , since this is an 
essential property of such materials [2] . We follow Ref. [6] 
and consider the following frequency dependencies of the 
effective left-handed medium characteristics, 



e 2 (w) = 1 - 



/U 2 (w) = 1 - 



Fu 2 



,2 ' 



(7) 



where the parameters correspond to the experimental 
data of Ref. [6]: uj p = 10 GHz, to = 4 GHz, and 
F = 0.56. The region of simultaneously negative permit- 
tivity and permeability in this case ranges from 4 GHz 
to 6 GHz. Dispersion curves for the first three guided 
modes in a slab waveguide with the thickness parameter 
L = 2 cm are shown in Fig. 2, where dashed curves cor- 
respond to "fast" modes, and solid - to "slow" modes. 
We find that the fundamental "slow" mode exists only at 



higher frequencies, whereas the second-order fast mode 
appears at lower frequencies. Both modes can have ei- 
ther positive or negative group- velocity dispersion in dif- 
ferent parameter regions. Properties of the first-order 
antisymmetric mode are different. The type of this mode 
seamlessly changes from "fast" to "slow" as the frequency 
grows. This transition occurs when the condition fc 2 = 
is satisfied, which is a boundary separating the two types 
of modes, as shown in Fig. 2 by a dotted line. The high 
"fast" modes exists at the frequencies close to the reso- 
nance at to = AGHz. 

In the left-handed materials, the electromagnetic waves 
are backward, since the energy flux and wave vector have 
opposite directions [2], whereas these vectors are paral- 
lel in conventional (right-handed) homogeneous materi- 
als. The energy flux is characterized by the Pointing vec- 
tor averaged over the period T — 2tt/lo and defined as 
S = (c/Sw) Re [E x H*]. A monochromatic guided mode 
has, by definition, a stationary transverse profile, and 
the averaged energy flux is directed along the waveguide 
only. It follows from the Maxwell's equations and Eq. (2) 
that the z-component of the energy flux is, 



8irujp(x) °' 



(8) 



The total power flux through the waveguide core and 
cladding can be found as P 2 = J^ L S z dx and Pi = 

2 S z dx, respectively. We find that the energy flux 
distribution for the waves guided along the left-handed 
slab is rather unusual. Indeed, the energy flux inside 
the slab (with p < 0) is opposite to that in the sur- 
rounding medium (with p > 0). This occurs because 
the normalized wave vector component along the waveg- 
uide (h) is fixed in a guided mode according to Eq. (2). 
An important information about the guided modes can 
be extracted from the study of the normalized energy 
flux P = (Pi + P 2 )/(|Pi| + |P 2 |). This parameter is 
bounded, |P| < 1, P — > 1 when the mode is weakly 
localized (|Pi| ^> IP2I), whereas P < for modes which 
are highly confined inside the left-handed slab. 

We have performed a detailed analysis of the slow 
guided modes and identified four distinct cases. 

(i) e 2 /i 2 > ei/ii, /i 2 > H\. Only odd mode exists below 
the threshold, p < p\/ p 2 - The corresponding critical 
value of the slab thickness L below which the odd mode 
exists is found as 



j _ c 



Ml 



LO p 2y Je 2 p 2 - £1^1 



(9) 



The energy flux P is positive for all values of L. The 
modes are forward propagating, i.e. the total energy flux 
along the waveguide is co-directed with the wavevector. 

(ii) e 2 p 2 > eipi, p 2 < p\. Even mode exists at all 
values of p, however odd modes can appear only when 
a threshold parameter value is exceeded, p > p 2 /p 2 . 
Accordingly, the critical value (9) determines the lower 



4 



X 



.AAAAAAAAAWMWrrjAAAAAAAi 
iUAAAAAAHTTTHTTTtAAAAlAii 

**mJZI! UIW- - - 



FIG. 4: Structure of the Pointing vector field in a localized 
surface wave propagating along a left-handed slab. 

boundary of the existence region for odd modes. The to- 
tal energy flux is negative for all L, and the modes are 
backward. The energy is mostly localized inside the slab. 

(hi) £2/i2 < eiMi; M2 > Ml- Both odd and even modes 
exist at all values of p and L, and the modes are forward. 

(iv) e2/i2 < eiMi: < Mi- Only even modes exist 
below the threshold value of p which can be found nu- 
merically using Eq. (5). Characteristic dependences of 
the wavenumber and normalized power on the slab width 
is shown in Figs. 3(a,b). At any slab thickness below a 
critical value, two modes always co-exist. One of the 
modes is forward and weakly localized, but the other one 
is backward and more confined. When the slab width ap- 
proaches the critical value, the branches corresponding to 
different modes merge, and the energy flux vanishes. In 
this special case, the energy fluxes inside and outside the 
slab exactly compensate each other. 

Since the energy fluxes are oppositely directed inside 
the guided modes, it might initially seem that such waves 
can only be sustained by two continuously operating 



emitters positioned at the opposite ends of the waveg- 
uide. Therefore, it is of the fundamental importance to 
understand whether wave packets of unite temporal and 
spatial extension can exist in left-handed waveguides. We 
calculate the Pointing vector averaged over the period 
of the carrier frequency, and present the characteristic 
structure of the energy flow in Fig. 4. Due to the unique 
double-vortex structure of the energy flow, most of the 
energy remains inside the wavepacket, and it does not 
disintegrate. The group velocity is proportional to the 
total energy flux P, and it can therefore be made very 
small or even zero by a proper choice of the waveguide 
parameters as demonstrated above. On the other hand, 
the group-velocity dispersion, which determines the rate 
of pulse broadening, can also be controlled. This flexibil- 
ity seems very promising for potential applications. 

Finally, we note that recent numerical simulations 
demonstrated that the phenomenon of the negative re- 
fraction, similar to that found for the left-handed meta- 
materials, can be observed in photonic crystals [8, 11]. 
Although in this case the wavelength is of the same or- 
der as the period of the dielectric structure (and, there- 
fore, a simple analysis in terms of the effective medium 
approximation is not strictly justified), we expect that 
similar mechanisms of wave localization will remain gen- 
erally valid. 

In conclusion, we have described, for the first time to 
our knowledge, linear guided waves in left-handed slab 
waveguides. We have demonstrated a number of exotic 
properties of such waves, including the absence of funda- 
mental modes and the sign-varying energy flux, and we 
have predicted the existence of the fundamentally novel 
classes of guided waves with a vortex-type internal struc- 
ture. 

We thank C. T. Chan and C. M. Soukoulis for use- 
ful discussions, and to D. E. Edmundson for assistance 
with Fig. 4. The work was partially supported by the 
Australian Research Council. 



* URL: http://www.rsphysse.anu.edu.au/nonlinear 
[1] R.A. Shelby, DR. Smith, and S. Shultz, Science 292, 77 
(2001). 

[2] V.G. Veselago, Usp. Fiz. Nauk 92, 517 (1967) [Sov. Phys. 
Usp. 10, 509 (1968)]. 

[3] M.C.K. Wiltshire, J.B. Pendry, I.R. Young, D.J. Lark- 
man, D.J. Gilderdale, and J.V. Hajnal, Science 291, 849 
(2001). 

[4] J.B. Pendry, A.J. Holden, W.J. Stewart, and I. Youngs, 

Phys. Rev. Lett. 76, 4773 (1996). 
[5] P. Markos and CM. Soukoulis, Phys. Rev. B 65, 033401 

(2001); D.R. Smith, S. Schultz, P. Markos, and CM. 

Soukoulis, Phys. Rev. B 65, 195104 (2002). 
[6] D.R. Smith, W. Padilla, D.C Vier, S.C Nemat-Nasser, 

and S. Shultz, Phys. Rev. Lett. 84, 4184 (2000). 



[7] J.B. Pendry, Phys. Rev. Lett. 85, 3966 (2000). 

[8] C. Luo, S.G. Johnson, JD. Joannopoulos, and J.B. 

Pendry, Phys. Rev. B 65, 201104 (2002). 
[9] R. Ruppin, Phys. Lett. A 277, 61 (2000); J. Phys. 

Condens. Matter 13, 1811 (2001); F.D.M. Haldane, 

arXiv:cond-mat:0206420 (2002). 
[10] M. B. Vinogradova, O. V. Rudenko, and A. P. Sukho- 

rukov, The Theory of Waves (Nauka, Moscow, 1990) (in 

Russian) . 

[11] H. Kosaka, T. Kawashima, A. Tomita, M. Notomi, T. 
Tamamura, T. Sato, and S. Kawakami, Phys. Rev. B 58, 
10096 (1998); M. Natomi, Opt. and Quant. Electron. 34, 
133 (2002). 



